Granular transport in a horizontally vibrated sawtooth channel 



Shahin Mobarakabadi/ Ehsan Ncdaacc Oskocc,^'El Matthias Schrotcr,^ and Mchdi Habibi^ 

'institute for Advanced Studies in Basic Sciences, Gava Zang, Zanjan 45195-159, Iran 
^Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Gottingen, Germany 

(Dated: February 25, 2013) 

We present a new mode of transport of spherical particles in a horizontally vibrated channel with 
sawtooth shaped side walls. The underlying driving mechanism is based on an interplay of directional 
energy injection transformed by the sidewall collisions and density dependent interparticle collisions. 
Experiments and matching numerics show that the average particle velocity reaches a maximum at 
60 % of the maximal filling density. Introducing a spatial phase shift between the channel boundaries 
increases the transport velocity by an order of magnitude. 

PACS numbers: 45.70.-n, 45.70.Mg, 83.80.Fg, 87.15. hj 



INTRODUCTION 



Brownian motors extract useful work from a noisy en- 
vironment by means of a broken spatial symmetry [H 0] ■ 
The concept of Brownian motors is of great importance 
in cell biology and nanotechnology, however their un- 
derlying principle is not limited to thermal noise and 
can therefore be also implemented in macroscopic se- 
tups. One example for such an athermal, macroscopic, 
and noisy environment arc granular gases 

Due to the dissipative collisions between particles, 
granular gases require constant external driving which is 
in most cases provided by shaking the container. Work 
can then e.g. be extracted by means of a rotational 
ratchet, where the symmetry is broken by different coat- 
ings on the two sides of eacli vane [1] . Or probe particles 
with an asymmetric shape can be set into translational 

Q or rotational Q motion. 

Another class of granular Brownian motors converts 
the random motion of the particles into a directed flow. 
This can either be implemented by breaking the sym- 
metry of the driving [3, or by breakin g th e spatial 
symmetry of the container boundaries [T^-Ufll- All se- 
tups in the latter group include a sawtooth shaped base 
plate which is shaken vertically to drive the granular gas. 
Their phenomenology includes height-dependent flow di- 
rections [lol - [l^ , segregation of binary mixtures [l^, [13] , 
and rotational motion of the circular base (l6j . 

In this paper we decouple the direction of driving from 
gravity: our system is a horizontally vibrated channel 
where sawtooth shaped side walls break the symme- 
try. We observe unidirectional transport with a non- 
monotone density dependence. In section 2 we discuss 
the experimental setup and results. In section 3 we in- 
troduce and validate our numerical simulations which are 
then used in section 4 to gain insight into the transport 
mechanisms. 
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FIG. 1: (Color online) Sawtooth channel. Top view (a), and 
dimensions (b) of the experimental setup and visualization of 
the setup used in the simulations (c). 



II. EXPERIMENTAL RESULTS 

The experimental setup consists of two sawtooth corri- 
dors which are connected by two small pentagonal cells at 
each end of the corridors, as shown in Fig. (Ta] Each cor- 
ridor consists of five triangular cells; both floor and walls 
are made from plexiglass. The channels are mounted 
horizontally and shaken in X direction (c.f. fig. Ilbp us- 
ing an off center pulley which is driven by an AC motor. 
The oscillations amplitude is 7.5 ± 0.5 mm, the driving 
frequency 5 Hz. 

The channel is filled with steel spheres of radius 6 mm. 
The maximum number of spheres per triangular cell is 
21, this value defines the maximal area fraction <f>c- Ten 
percent of the beads arc marked with a different color. 
The average speed of the beads is determined by measur- 
ing the rate at which tracers enter a given cell; the sample 
time is chosen so that the average tracer has completed 
at least two full cycles. 

Figure [5] shows the average drift speed as a function of 
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FIG. 2: Average drift speed as a function of the reseated area 
fraction, measured botli in experiment and simulation. Ver- 
tical dashed lines indicate the transition between the dilute 
(left), cooperative (middle) and jamming regime (right). The 
three example error bars are standard deviations over 10 cy- 
cles. The simulation time corresponds to the average particle 
performing four cycles. 

the renormalized area fraction (4)/(j)c) where (j) is the num- 
ber of spheres per cell. Three different regimes can be dis- 
tinguished. For low area fractions (0.005 < (j)/4>c < 0.08) 
the particles move in negative Y direction (i.e. downward 
in fig. Ilb[) with a velocity that is independent of (j). As 
here are on average less than 2 beads in a cell, particle- 
particle collisions arc rare. We refer to this range as the 
dilute regime. 

For intermediate area fractions (0.08 < (f>/(t)c < 0.6) 
the average speed increase monotonously. Fitting a 
line to the experimental data in figure [2] shows that 
i^y) ^ {4>l4'cf''^ in the intermediate regime. This effect 
delineates the cooperative regime. 

For even higher densities {4)/4>c > 0.6) the drift speed 
decreases sharply. As the system approaches the fluid- 
solid transition density, we refer to this as the jamming 
regime. 



III. SIMULATION METHOD 

To elucidate the transport mechanism we have per- 
formed molecular dynamics simulations using a simpli- 
fied version of the contact force laws proposed by Silbert 
et al. . The normal and tangential forces between two 
grains i and j, F„.^ and Ft^^., and are given by 

S' 

= /(-^)(fc«^i,jnij - 7nme//V„, J (1) 

and 

Fti,3 = ~ niin(7tTOe//Wt,^ , /isi^„i,j )vt,^^ , (2) 
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TABLE I: The numerical values of the contact model coeffi- 
cients. The left column refers to particle-particle interactions, 
the right column to particle-wall. All values are given in the 
CGS system of units. 



Here v„ . ^ and v^^ . are the normal and tangential com- 
ponents of the relative velocity between the grains (eval- 
uated at the contact point) . . ^ is the unit vector of tan- 
gential velocity, n; j- the unit vector corresponding to the 
distance r^j between the two grain centers, d is the par- 
ticle diameter, rueff the effective mass miirij / (rrii + mj), 
and 6ij ~ d — ri^ is the normal compression at the con- 
tact. f{x) equals ^/x which models Hertzian contacts. 
The density of the grains is set to 7.87 g/cm'^. The nu- 
merical values of the elastic constant kn , the viscoclastic 
constants 7„ and 74 , and the coefficient of friction fi are 
listed in table HI 

We do neglect the the viscoclastic tangential interac- 
tion between grains by setting 74 = 0, this assures that 
particles do not experience any vertical forces and remain 
therefore in contact with the bottom surface of the chan- 
nel. The agreement shown below between our simulations 
and experiments justifies this decision by hindsight. 

The same equations are used to model the interaction 
between grains and the channel walls, the coefficients 
have an added subscript w. The mass of the container is 
set to be infinity, consequentially mefj = m. 

For the evaluation of the tangential force f exerted by 
the channel bottom on a particle we follow a method pro- 
posed by Kondic [lH which distinguishes between sliding 
and rolling contacts. In a first step the contact is assumed 
to be rolling with f = ma. Taking the rotation into ac- 
count, the acceleration a of the sphere can be computed 
as a = 2/7as with a^. being the acceleration of the sur- 
face. If the no-sliding condition (|f | < fmax = ^-sfng) is 
satisfied, this is the final result. 

Otherwise, the sliding contact has to be evaluated us- 
ing the dynamic friction coefficient ^dw'- 

f = -tJ-dw-mg^^ (3) 

where v^, is the velocity of the contact point relative to 
the substrate. 

The shaking of the container is modeled by a sinu- 
soidal excitation x{t) = xq cos (tut), uj is the oscillation 
frequency and xq is its amplitude. The code is written in 
C-|--|- using a 5th order predictor-corrector algorithm for 
numerical integration of the equations of motion. The 
time step increment (At) is set to 10~^ in all simula- 
tions; the simulation time is chosen so that the average 
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FIG. 3: Modified boundary geometries. Channel a) interpo- 
lates between a sawtooth and a triangular shape by increasing 
the angle a. /3 is kept constant at 72 degree as in figure [T] 
Channel b) introduces a relative displacement d of the two 
sidewalls. 

particle has performed four full cycles which requires at 
maximum 3.6 x 10^ time steps. 

Figure [5] shows that the results of simulation are in 
good agreement with the experiment. Therefore we will 
use our simulations to gain insight into the driving mech- 
anism in the three regimes. 

IV. TRANSPORT MECHANISM 

Understanding the transport mechanism depicted in 
figure [2] requires answers to the following three questions: 

(A) Why are the particles moving unidirectionally at all? 
This is best discussed in the dilute regime. 

(B) Why does the drift velocity increase in the coopera- 
tive regime? 

(C) Why does the drift velocity go down at even higher 
densities? 



A. Particles moving in the dilute regime 

The dilute regime is specifically apt to study the un- 
derlying transport mechanism as it allows us to ignore 
inter-particle collisions. First, it should be pointed out 
that dissipative collisions are not a necessary condition 
to have a finite drift velocity; the mechanism is therefore 
different from fluxes seen in granular gases going through 
beveled pores [l^. This has been tested by short simu- 
lations where 7„ and jnw have been set to zero. 



FIG. 4: The average drift speed depends on the asymmetry 
of the channel. Data are measured in setup |3a] at a shaking 
frequency of 6.4 Hz. 

Secondly, like in other ratchet systems the symmetry 
breaking of the boundary is a necessary condition to ob- 
tain a finite (vy). If we increase in our simulations the 
angle a shown in Figure I3a[ we change the channel pro- 
file from a sawtooth shape to a triangular one. This goes 
together with a reduction of the average drift velocity 
towards zero as shown in figure 51 

We therefore conjecture the following driving mech- 
anism: spheres gain kinetic energy in x-direction due 
to the frictional contacts with the horizontally vibrat- 
ing base. This kinetic energy is then converted by the 
sidewall collisions into a directed motion in negative y- 
direction. To substantiate this conjecture we have com- 
puted the distribution P of particle direction with respect 
to the Y-axis 7, measured after the particles collided with 
the sidewall sidewall (c.f. Fig. Ilb[) . As it can be seen in 
figure[5]the majority of post coUisional velocities point in 
negative y-direction which corresponds to 7 values larger 
than 90 degree. For (p/cj)^ = 0.047 more than 60 % of 
post collision vectors fall into this range. 



B. Cooperative regime 

For values of (p/ipc > 0.08, the average drift velocity 
increases proportional to ((/)/0c)°'® as shown in figure 
[21 However, the post sidewall collision angle distribution 
P{j) even becomes more balanced, at 4>/4>c = 0.476 only 
58 % of the vectors point in negative y-direction. There- 
fore the mechanism behind the increase of (uy) has to be 
something different. We hypothesize that the increased 
density keeps the particles longer in close proximity to 
the narrow passage between cells. 

There are two indications for the importance of the 
narrow passage area. The first one can be seen in figure 
iniwhich shows that (vy) has (at all densities) a local max- 
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FIG. 5: Angular distribution of particle velocities after colli- 
sions with the sawtooth side walls. The angle 7 is defined in 
figure [Tbl Values of 7 between 90 and 198 degrees make the 
particles go into the direction of the average drift. 



FIG. 7: Shifting the two sidewall by a distance d/ L as dis- 
played in figure [3bl increases the average drift velocity by an 
order of magnitude and removes the transition between the 
dilute and the cooperative regime. 



imum when the shaking amplitude xq is approximately 
half of the width of the passage. It is also interesting 
to note that only in the in the cooperative regime the 
average drift speed increases systematically with the am- 
plitude Xq. 

The second sign that improved overcoming of the nar- 
row passage is the mechanism behind the cooperative 
regime can be found when studying a channel where the 
two sidewalk are shifted by a length d against each other 
as in figure I3bl Figure [7] clearly demonstrates that the 
distinction between dilute and cooperative regime van- 
ishes when the ratio d/L goes to 0.5 i.e. the narrow pas- 



sages arc taken out of the system. The accompanying 
increase of the average drift speed by one order of mag- 
nitude is also significant from the perspective of potential 
applications. 



C. Jamming regime 

The strong decrease in drift velocity for (p/cjic > 0.6 can 
be traced back to two mechanisms. First the effective 
redirection of the kinetic energy at the sidewalls stops. 
As can be seen in figure Oat (f>/<f>c = 0.857 only 49 % of 
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FIG. 6: Average drift speed versus shaking amplitude for 
three different values of cj>/4>c- The drift amplitude increases 
only in the cooperative regime with the driving strength. 
Measured in the standard setup displayed in figure [Tc] 



FIG. 8: The average collision duration doubles in the jamming 
regime, indicating a strong decrease of kinetic energy. Data 
are measured in the standard setup (figure [Tc)) at a shaking 
frequency of 6.4 Hz. 
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the post coUisional velocities have 7 values larger than 
90 degree. 

Secondly, figure [8] shows that the average collision du- 
ration doubles in the jamming regime. For Hertzian par- 
ticles the collision duration is proportional to v~^'^ j20| . 
Therefore the average kinetic energy of the particles be- 
comes roughly a factor 1000 smaller for (l)/4'c increasing 
from 0.6 to 0.9. At the same time the number of (dissipa- 
tive) collisions per time only doubles. We conclude that 
the energy injection through the shaking base plate be- 
comes much less efficient, probably due to the formation 
of short-lived but system-spanning force chains. 

V. CONCLUSIONS 

We have demonstrated experimentally and numerically 
a new method of granular transport in a vertically vi- 
brated channel with sawtooth shaped boundaries. The 
average drift speed depends on the filling fraction, in- 



creasing for intermediate densities and decreasing again 
when approaching jamming. The underlying driving 
mechanism was found to be the redirection of the in- 
jected kinetic energy in the particle sidewall collisions. 
The reason for the drift speed increase at intermediate 
densities seems to be the improved fiux through the nar- 
row passages between the tips of the sawteeth; remov- 
ing them by introducing a phase shift between the two 
boundaries resulted in a more than tenfold increase in 
drift velocity. At high densities the average contact du- 
ration increases which indicates that the energy injection 
through the shaking bottom plate becomes significantly 
less efficient. However, more work will be needed for 
establishing a microscopic theory of this transport mech- 
anism. 
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